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Abstract. We extend methods of Greenberg and the author to compute 
in the cohomology of a Shimura curve defined over a totally real field 
with arbitrary class number. Via the Jacquet-Langlands correspondence, 
we thereby compute systems of Hecke eigenvalues associated to Hilbert 
modular forms of arbitrary level over a totally real field of odd degree. 
We conclude with two examples which illustrate the effectiveness of our 
algorithms. 



The development and implementation of algorithms to compute with au- 
tomorphic forms has emerged as a major topic in explicit arithmetic geometry. 
The first such computations were carried out for elliptic modular forms, and now 
very large and useful databases of such forms exist |2I13I14) . Recently, effective 
algorithms to compute with Hilbert modular forms over a totally real field F 
have been advanced. The first such method is due to Dembele |4l5j , who worked 
initially under the assumption that F has even degree n — [F : Q] and strict 
class number 1. Exploiting the Jacquet-Langlands correspondence, systems of 
Hecke eigenvalues can be identified inside spaces of automorphic forms on i?^, 
where B is the quaternion algebra over F ramified precisely at the infinite places 
of F — whence the assumption that n is even. Dembele then provides a compu- 
tationally efficient theory of Brandt matrices associated to B. This method was 
later extended (in a nontrivial way) to fields F of arbitrary class number by 
Dembele and Donnelly [5]. 

When the degree n is odd, a different algorithm has been proposed by Green- 
berg and the author [5] , again under the assumption that F has strict class num- 
ber 1. This method instead locates systems of Hecke eigenvalues in the (degree 
one) cohomology of a Shimura curve, now associated to the quaternion algebra 
B ramified at all but one real place and no finite place. This method uses in a 
critical way the computation of a fundamental domain and a reduction theory 
for the associated quaternionic unit group [TC] ; see Section 1 for an overview. In 
this article, we extend this method to the case where F has arbitrary (strict) 
class number. Our main result is as follows; we refer the reader to Sections 1 
and 2 for precise definitions and notation. 



Theorem 1. There exists an (explicit) algorithm which, given a totally real field 
F of degree n = [_F : Q] , a quaternion algebra B over F ramified at all but one 
real place, an ideal 71 of F coprime to the discriminant 2) of B, and a weight 
k € (2Z>o)", computes the system of eigenvalues for the Hecke operators Tp with 
p I J)*Jt and the Atkin-Lehner involutions Wpc with p*^ \\ DVl acting on the space 
of quatcrnionic modular forms {'^1) of weight k and level for B. 

In other words, there exists an exphcit finite procedure which takes as input 
the field F, its ring of integers Zp, a quaternion algebra B over F, an ideal 
0^ C Zi?, and the vector k encoded in bits (each in the usual way), and outputs 
a finite set of number fields Ef C Q and sequences (a/(p))p encoding the Hecke 
eigenvalues for each cusp form constituent / in 5^(5^1), with a/(p) € Ef. 

From the Jacquet-Langlands correspondence, applying the above theorem to 
the special case where D = (1) (and hence n = [F : Q] is odd), we have the 
following corollary. 

Corollary 2. There exists an algorithm which, given a totally real field F of odd 
degree n — [F : Q], an ideal 9T of F, and a weight k € (2Z>o)", computes the 
system of eigenvalues for the Hecke operators Tp and Atkin-Lehner involutions 
Wpc acting on the space of Hilbert modular cusp forms 5^(91) of weight k and 
level 0^. 

This corollary is not stated in its strongest form: in fact, our methods overlap 
with the methods of Dembele and his coauthors whenever there is a prime p 
which exactly divides the level; see Remark [5] for more detail. Combining these 
methods, Donnelly and the author [7^ are systematically enumerating tables of 
Hilbert modular forms, and the details of these computations (including the 
dependence on the weight, level, and class number, as well as a comparison of 
the runtime complexity of the steps involved) will be reported there [7], after 
further careful optimization. 

A third technique to compute with automorphic forms, including Hilbert 
modular forms, has been advanced by Gunnells and Yasaki |^. They instead use 
the theory of Voronoi reduction and sharbly complexes; their work is independent 
of either of the above approaches. 

This article is organized as follows. In Section 1, we give an overview of the 
basic algorithm of Greenberg and the author which works over fields F with 
strict class number 1. In Section 2, using an adelic language we address the 
complications which arise over fields of arbitrary class number, and in Section 3 
we make this theory concrete and provide the explicit algorithms announced in 
Theorem [T] Finally, in Section 4, we consider two examples, one in detail; our 
computations are performed in the computer system Magma [I]. 

The author would like to thank Steve Donnelly and Matthew Greenberg for 
helpful discussions as well as the referees for their comments. The author was 
supported by NSF Grant No. DMS-0901971. 



1 An overview of the algorithm for strict class number 1 



In this section, we introduce the basic algorithm of Greenberg and the author 
[8] with a view to extending its scope to base fields of arbitrary class number; 
for further reading, see the references contained therein. 

Let F be a totally real field of degree n — [F : Q] with ring of integers lip- Let 
Fj^ be the group of totally positive elements of F and let Zp = n F^ . Let 
-B be a quaternion algebra over F of discriminant D . Suppose that B is split at a 
unique real place vi , corresponding to an embedding ioo : B S (8) M = M2 (M) , 
and ramified at the other real places V2, ■ ■ ■ ,Vn- Let 0(1) C B be a maximal 
order and let 

0{l)l - {7 e 0{iy^ : ^i(nrd(7)) > 0} - {7 G 0{1) : nrd(7) G 

denote the group of units of 0(1) with totally positive reduced norm. Let 

r(l) = /,oo(0(l)^/Z^) C PGL2(M)+, 

so that r{l) acts on the upper half-plane H ^ {z G C : Im(z) > 0} by linear 
fractional transformations. Let *Tt C Zj? be an ideal coprime to D, let O = Co(^) 
be an Eichler order of level m, and let F = Foim) = too(Co(^)^ /Z^). 

Let k = {ki, . . . ,kn) € (2Z>o)" be a weight vector; for example, the case 
k = (2, ... ,2) of parallel weight 2 is of significant interest. Let 5*^(0^) denote 
the finite-dimensional C- vector space of quaternionic modular forms of weight k 
and level 91 for B. Roughly speaking, a form / €E S^{Ul) is an analytic function 
f : "H ^ Wk(C) which is invariant under the weight k action by the group 7 € -T, 
where Wfe(C) is an exphcit right B^-module [H (2.4)] and Wk{C) = C when k 
is parallel weight 2. The space S^{Ul) comes equipped with the action of Hecke 
operators Tp for primes p \ 2)91 and Atkin-Lehner involutions Wpi for prime 
powers II DDI. 

The Jacquet-Langlands correspondence 'W, Theorem 2.9] (see Hida [lOl Propo- 
sition 2.12]) gives an isomorphism of Hecke modules 

5f (91) A 5fe(D91)®-"°"', 

where S'fc(J'Ol)®"'^''™ denotes the space of Hilbert modular cusp forms of weight k 
and level 2)91 which are new at all primes dividing S. Therefore, as Hecke mod- 
ules one can compute equivalently with Hilbert cusp forms or with quaternionic 
modular forms. 

We compute with the Hecke module 5*^(91) by identifying it as a subspace 
in the degree one cohomology of -r(l), as follows. Let T4(C) be the subspace 
of the algebra <C[xi,yi, . . . ,a:„,?/„] consisting of those polynomials q which are 
homogeneous in (xi,yi) of degree Wi = ki — 2. Then Vfc(C) has a right action of 
the group B" given by 

q'^{xi,yi,...,Xn,y„) = ^ JJ(det 7,)""''/^ j q{{xi 2/1)71, (a;„ ?/„)7n) (1) 



for 7 e B^, where ~ denotes the standard hwolution (conjugation) on B and 
7i — Vi{'j) £ M2(C). By the theorem of Eichler and Shimura [51 Theorem 3.8], 
we have an isomorphism of Hecke modules 

Ai/i(r,Ffe(C))+ 

where the group cohomology denotes the (finite-dimensional) C- vector space 
of crossed homomorphisms f : F ^ Vfe(C) modulo coboundaries and + denotes 
the -|-l-eigenspace for complex conjugation. By Shapiro's lemma [51 §6], we then 
have a further identification 

s^m ^ H'{r,Vk{C))^ = H\r{i),vic))+, (2) 

where ^(C) = Coind^^^ 14(C). 

In the isomorphism dSj), the Hecke operators act as follows. Let p be a prime 
of Zi? with p I and let Fp denote the residue class field of p. Since F has 
strict class number 1, by strong approximation [15[ Theoreme III. 4. 3] there exists 
77 € C such that nrdTr is a totally positive generator for p. It follows that there 
are elements 7^ G O^, indexed by a S P^(Fp), such that 

0^7tO^= □ Olaa (3) 

aGPi(Fp) 

where aa = irja- 

Let / : r{l) — > F(C) be a crossed homomorphism, and let 7 G -^(1). The 
decomposition ^ extends to 0(1) as 

Oil)lnOil)l= □ 0(1)>.. 

Thus, there are elements 5a G 0(1)+ for a G IP^(Fp) and a unique permutation 
7* of pi(Fp) such that 

aa7 = SaO^-a (4) 

for all a. We then define / 1 Tp : r{l) V{C) by 

aePi(Fp) 

The space 5*^(91) similarly admits an action of Atkin-Lehner operators Wpe for 
primes p*^ || DOT. 

From this description, we see that the Hecke module H^{r{l),V{C))~^ is 
amenable to explicit computation. First, we compute a finite presentation for 
r{l) with a minimal set of generators G and a solution to the word problem for 
the computed presentation using an algorithm of the author [TB]. Given such a 
set of generators and relations, one can explicitly find a basis for the C-vector 
space H\r{l),ViC)) [i §5]. 



Wc then compute the action of the Hecke operator Tp on H^{r{l),V{C)). 
We first compute a sphtting Lp : O ^ M2{'Zp^p). The elements aa in ^ are then 
generators with totally positive reduced norm of the left ideals 

Ia=OLp'(llj+Op (6) 

and are obtained by principalizing the ideals la', here again we use strong ap- 
proximation and the hypothesis that F has strict class number 1. Then for each 
a £ P^(Fp) and each 7 e G, we compute the permutation 7* [8l Algorithm 
5.8] and the element Sa = aaja^^^ e r{l) as in (|4]). Using the solution to the 
word problem, we then write Sa as a word in the generators G for r{l), and 
then for a basis of crossed homomorphisms / we compute f \Tp by computing 
(/ I ^p)(7) ^ V{C) for each 7 G G as in ([5]). In a similar way, we compute the 
action of complex conjugation and the Atkin-Lehner involutions. We then de- 
compose the space H^{r, V{C)) under the action of these operators into Hecke 
irreducible subspaces, and from this we compute the systems of Hecke eigenval- 
ues using linear algebra. 

2 The indefinite method with arbitrary class number 

In this section, we show how to extend the method introduced in the previous 
section to the case where F has arbitrary class number [8l Remark 3.11]. We 
refer the reader to Hida [11] for further background. 

2.1 Setup 

We carry over the notation from Section 1. Recall that O — Oq{'^) is an Eichler 
order of level *Jt in the maximal order 0(1) C B. 

Let = {z € C : Im(z) ^ 0} = C \ M be the union of the upper and 
lower half-planes. Then via too, the group i?^ acts on by linear fractional 
transformations. 

In this generality, we find it most elucidating to employ adelic notation. Let 
Z — ^im Z/nZ and let ^ denote tensor with Z over Z. Consider the double coset 

x(c) = B^\(H± X i?v<5''), 

where B^ acts on B^/O^ by left multiplication via the diagonal embedding. 
Then X{C) has the structure of a complex analytic space |;31 which fails to be 
compact if and only if i? = M2 (Q) , corresponding to the classical case of elliptic 
modular forms — higher class number issues do not arise in this case, so from 
now we assume that i? is a division ring. 

We again write S^{m) for the finite-dimensional C- vector space of quater- 
nionic modular forms of weight k and level here, again roughly speaking, 
a quaternionic modular form of weight k E (2Z>o)" and level 91 for B is an 
analytic function 

f -.n^x B^'/d' Wk{C) 
which is invariant under the weight k action of B^ , with Wk{C) as in Section 1. 



2.2 Decomposing the double coset space 

By Eichler's theorem of norms, we have md(-B^) = F^_^^ where 

Fj^^j = {a e F"" : v^{a) > for i = 2, . . . , n} 

is the subgroup of elements of F which are positive at all real places which are 
ramified in B. In particular, /B^ = Z/2Z, where 

= {7 e : t.i(nrd(7)) > 0} = {7 G B : nrd(7) € F^}. 

The group B^ acts on the upper half-plane H, therefore we may identify 

X{C) = Bl\{H X B^'/d''). 
Now we have a natural (continuous) projection map 

x{c) B^\i?v^^ 

and by strong approximation [ 151 Theoreme III. 4. 3] the reduced norm gives a 
bijection 

nrd:B^\BVO'' ^^^\^V% = C1+Zf, (7) 

where CI'*' Zp denotes the strict class group of Z^?, i.e. the ray class group of Zp- 
with modulus equal to the product of all real (infinite) places of F. 

The space ^(C) is therefore the disjoint union of Riemann surfaces indexed 
by CY^Zp, which we identify explicitly as follows. Let the ideals b G Zp form 
a set of representatives for Cl^ Z^, and let 6 G Z^? be such that bZp (iZp = b. 
For expositional simplicity, choose b = Zp and /3 — 1 for the representatives of 
the trivial class. By strong approximation ([7]), there exists (3 E B^' such that 
nrd(/3) = b. Therefore 

xic) = \jB^inxpd''). (8) 

lb] 

We have a map 

where = (30/3-^ n B and =0^nBl, so that = O. 

For each P, let F^ = l^{0^ JZ'^) C PGL2(K)+. Then the Eichler-Shimura 
isomorphism on each component in (|5]) gives an identification of Hecke modules 

S^im)^^H'iF^,V,iC))+, (9) 

where + denotes the -t-l-eigenspace for complex conjugation. For each /?, let 
0(1)3 — I30{l)f3~^ n B he the maximal order containing the Eichler order O3, 



and let r(l)^ = i^{0{l)^ JZ^'p). Further, let ^^(C) = Coind^^^^^^^^ Vu{C). Then 
Shapiro's lemma applied to each summand in ([9]) gives 

s^m ^ 0ffi(r(i)^,F^(c))+. (10) 

2.3 Hecke operators 

In the description (fTO|) . the Hecke operators Tp act on 0^ _ff ^(/^(l)^, Vg(C)) in 
the following way. Let p be a prime ideal of Zp- with p \ D'Vl, and let p g "Lp be 
such that pl^F n Zi? = p. We consider the /3'-summand in PH)) . corresponding 
to the ideal class [b']. Let / : -^(1)^/ ~^ Vp{'C) be a crossed homomorphism: we 

will then obtain a new crossed homomorphism f\Tp : ril)-^ — > V^^(C), where /3 
corresponds to the ideal class of [pb'] among the explicit choices made above. 

Let w G Op be such that nrd(TO) = p. Then there are elements 7a e O^, 
indexed by a G P^(Fp), such that 

a|ea|= u <5|s. (11) 

oepi(Fp) 

where = ^7a. 

Let 7 e rjj. Extending pT|) to 0(1) ^, we conclude that there exist unique 
elements 5a G C'(l)2 and a unique permutation 7* of P^(Fp) such that 

Sa7 = Saa^-a 

for a e pi(Fp). Thus we have 

where S'a^0'p-ya0'd-')-'- 

Recall that (3'0 has left order Op, and similarly 0/3~^ has right order Op. 

Therefore, we may consider the left O^j-ideal 

dp,d'dr^dpaa (12) 

noting that the left and right orders in each case match up, so the product is 
compatible. Next, recall that the elements /?', f3, w have reduced norms corre- 
sponding to the ideal classes [b'], [pb'], and [p], respectively. Thus the reduced 
norm of the left ideal has a trivial ideal class. Therefore, by strong approx- 
imation (applied now to left ideals of the order C^,), for each a E P^(Fp), there 
exist elements tt^ e Os, n such that 



d^,p'r'aanB^o^ya- 



Hence there exists a unique permutation 7* of P^(Fp) such that 

/ r/ / 

with 5a e O-, y The new crossed homomorphism / 1 Tp : ^ ''^(C) is then 
defined by the formula 

(/irp)(7)= E fi^'y^ 

for 7 e r^. 

2.4 Complex conjugation and Atkin-Lehner involutions 

We now define an operator Woo which acts by complex conjugation. Let Cl^^'' Zp 
denote the ray class group of Z^? with modulus equal to the real (infinite) places 
of F which are ramified in B. Then we have a natural map C1+ Zf -)■ Cl(+^ Zp; 
this map is an isomorphism if and only if there exists a unit u G Zp which 
satisfies vi{u) < and Vi{u) > for the other real places Vi [i = 2, . . . ,n) of 
F, otherwise the kernel of this map is isomorphic to Z/2Z. Let [m] € CV^ Zp 
generate the kernel of this map. 

Let / : ^(1)^/ — >■ V^, be a crossed homomorphism, and let (3 correspond to the 
ideal class [b'm~^]; we will define the complex conjugate crossed homomorphism 
(/ I VFoo) : r{l)^ V^{C). The left e'^,-ideal d^,P'd^-^nB has reduced norm 
corresponding to the ideal class [tn] G QX^Zp, so there exists a generator \i' e 
O-^, of this ideal such that wi(nrd(/x')) < but Wi(nrd(/i')) > for i = 2, . . . , n. 
Then given 7 G r(l)^, we define 

(/it^oo)(7) = /(MV-'r'- 

Finally, we define the Atkin-Lehner involutions W^c for p'^ || S)OT. Let p 
correspond to p G Z^? . Then there exists an element tt e which generates the 
unique two-sided ideal of O-^ of reduced norm generated by jf^ . The element tt 
normalizes O-^ and 7?^ e O^F^ . Let /3 correspond to the ideal class [pb']. Then 
as above, by strong approximation there exists an element /z' e O^, fl such 
that OpP'^Tf n B = Opii'. Given / : r(l)^, V^, we then define (/ \Wpe) : 

r(i)^^y^(C)by 

(/|W^p)(7) = /(MV-'r' 

for 7 e r(i)^. 

3 Algorithmic methods 

In this section, we take the adelic description of Section 2 and show how to 
compute with it explicitly, proving Theorem 1. 



Our algorithm takes as input a totally real field F of degree [F : Q] = n, 
a quaternion algebra B over F split at a unique real place, an ideal ^ d Ijf 
coprime to the discriminant S of S, a vector k G (2Z>o)", and a prime p | S^Tl, 
and outputs the matrix of the Hecke operator Tp acting on the space H = 

0^ _ff ^ (/^(l)^, V^{C))^ (in the notation of Section 2) with respect to some fixed 
basis which does not depend on p. From these matrices, one decomposes the space 
H into Hecke-irreducible subspaces by the techniques of basic linear algebra. 

Our algorithm follows the form given in the overview in Section 1, so we 
describe our algorithm in steps, with a description of each step along the way. 

Step 1 (Compute a splitting field): Let K ^ C he a Galois number field con- 
taining F which splits B: for example, we can take the normal closure of any 
quadratic field contained in B. Since all computations then occur inside K C C, 
we may work then with coefficient modules over K using exact arithmetic. (This 
step is only necessary if k is not parallel weight 2, for otherwise the action of 
factors through K ^ Q.) 

Step 2 (Compute ideal class representatives): Compute a set of representatives 
[b] for the strict class group Cl^ Zp with each b coprime to pSDI. (See Remark 
H] below.) 

Compute a maximal order 0(1) C B. For each representative ideal b, com- 
pute a right C'(l)-ideal J(, such that nrd(J(,) = & and let C(l)b be the left 
order of J^- (In the notation of Section 2, the right C'(l)-ideals Jt, represent the 
elements P, and 0{l)b = 0{1)^.) 

Step 3 (Compute presentations for the unit groups): Compute an embedding 
Lao '■ B ^ M2(M) corresponding to the unique split real place. 

For each b, compute a finite presentation for i"(l)b = ioo(0(l)^_^/Z^) con- 
sisting of a (minimal) set of generators Gf, and relations i?b together with a 
solution to the word problem for the computed presentation [TB]. (Note that 
the algorithm stated therein ^Gj Theorem 3.2] is easily extended from units of 
reduced norm 1 to totally positive units.) 

For efficiency, we start by computing such a presentation with generators G 
associated to the order 0(1) and then for each order 0(1)[, we begin with the 
elements in hand formed by short products of elements in G which happen to lie 
in O ( 1 ) f, (to aid in the search for units 16', Algorithm 3.2]; note that O ( 1 ) n O ( 1 ) t, 
is an Eichler order of level b in 0(l)(i). 

Step 4 (Compute splitting data): Compute a splitting 

: 0(1) ^ 0(1) Ip^m = M2{Zfm)- 

Note that since b is coprime to 91, we have 0(1) (E) = 0(l)b <E) I'P.'n for all 
b, so L<yi also gives rise to a splitting for each 0(1)[,. For each b, compute the 
Eichler order Ot, C 0(l)i, of level ^ with respect to 

Next, for each b, compute representatives for the left cosets of the group 
A = LooiO^ + f^F) inside r(l)i, 8, Algorithm 6.1]. Finally, identify 

V{K)i, = Coind;[;^^'^' Vk{K) 



as a K-vectoT space given by copies of Vk{K) indexed by these cosets, and 
compute the permutation action of the representatives of these cosets on this 
space. 

In practice, it is more efficient to identify the above coset representatives with 
elements of P^(Zi?/N) and thereby work directly with the coefficient module 

- K[F\ZF/m)] ® Vk{K). 

Step 5 (Compute a basis for cohomology): Identify the space of crossed homo- 
morphisms Z^{r{l)h, V{K)h) with its image under the inclusion 

Z\r,,V{K),)^ V{K), 
geG, 

consisting of those / G ©ggg^ V{K)ii which satisfy the relations f{r) = for r G 
Ri,. Compute the space of principal crossed homomorphisms B^{r{l)b, V{K)b) 
in a similar way, and thereby compute using linear algebra a A'-basis for the 
quotient H\r{l)b,V{K)b) = ZHr{l)b,V{K)b)/B\r{l)t„V{K)b) for each b. 
LetH = ®,H'{r{l),,V{K),). 

Step 6 (Compute representatives for left ideal classes): Compute a splitting tp : 
0(1) ^ M2{'Zf,p)- For each ideal b', perform the following steps. 

First, compute the ideal b with ideal class [b] — [pb']. Compute the left ideals 

indexed by the elements a — {x : y) E P^(Fp) and then compute the left Of- 
ideals 1'^ = Jb'Jbla- 

Compute totally positive generators tt^ G Of,' ^ for Ob'Tr^ = [H]. 

Now, for each 7 G Gfa, compute the permutation 7* of P^(Fp) [51 Algorithm 
5.8] and then the elements 5'^ — t:'^j'k'~^1^ for a G P"'^(Fp); write each such element 
(5^ as a word in G'^, and from the formula 

(/irp)(7)= E fi^'y^ 

with / in a basis for the b'-component of cohomology as in Step 5 compute the 
induced crossed homomorphism f \Tp in the b-component. 

Step 7 (Compute the blocks of the intermediate matrix): Assemble the matrix T 
with rows and columns indexed as in Step 5 with blocks in the (b, b') position 
given by the output of Step 6; this matrix describes the action of Tp on H. 

Step 8 (Decompose H into ±-eigenspaces for complex conjugation): Determine 
the representative ideal m (among the ideals b) which generates the kernel of 
the map C1+ Zp Cl^^^ Zp- 



For each ideal b', perform the foUowing steps. Compute the ideal b such 
that [b] = [b'm^^], and compute a generator /i' with Ob'fj.' = JwJb such that 
i)(nrd(/i')) < 0. For each 7 S Gj,, from the formula 

for / in a basis for the b'-component of cohomology as in Step 5 compute the 
induced crossed homomorphism / | Tp in the b-component. 

Assemble the matrix with blocks in the (b, b') position given by this output: 
this matrix describes the action of complex conjugation Woo on H . Compute a 
i^-basis for the +l-eigenspace 7?+ of H for Woo- Finally, compute the matrix 
T"*" giving the action of Tp restricted to and return T+. 

This completes the description of the algorithm. 

In a similar way, one computes the Atkin-Lehner involutions, replacing Step 6 
with the description given in Section 2.4, similar to the computation of complex 
conjugation in Step 8. 

Remark 3. Note that Steps 1 through 3 do not depend on the prime p nor the 
level ^ and Steps 4, 5, and 8 do not depend on the prime p, so these may be 
precomputed for use in tabulation. 

Remark 4- To arrange uniformly that the ideals b representing the classes in 
Cl^ Z p are coprime to the prime p in advance for many primes p , one has 
several options. One possibility is to choose suitable ideals b of large norm in 
advance. Another option is to make suitable modifications "on the fly": if p is 
not coprime to b, we simply choose a different ideal c coprime to p with [b] = [c], 
a new ideal Jc with nrd(Jc) = c, and compute an element v £ Ot, such that 
pOf)!^^^ = Oc- Conjugating by v where necessary, one can then transport the 
computations from one order to the other so no additional computations need 
to take place. 

4 Examples 

In this section, we compute with two examples to demonstrate the algorithm 
outlined in Section 3. Throughout, we use the computer system Magma [T]. 

Our first and most detailed example is concerned with the smallest totally 
real cubic field F with the property that the dimension of the space of Hilbert 
cusp forms of parallel weight 2 and level (1) is greater than zero and the strict 
class number of F is equal to 2. This field is given by F = Q{w) where w 
satisfies the equation f{w) — — llw — 11 = Q. The discriminant of F is equal 
to 2057 = ll^lT, and Zp = Z[w]. The roots of / in M are -2.602 . . . , -1.131 . . . , 
and 3.73 . . . , and we label the real places vi, U2, W3 of F into M according to this 
ordering. 

We define the sign of a € to be the triple sgn(a) = {sgn{vi{a)))^^i G 
The unit group of F is generated by the elements —1, w + 1 with sgn(u' + 1) = 
(1,-1,-1), and the totally positive unit —w^+2w + 12. 



We begin by finding a quaternion algebra B with T) ~ which is ramified at 
all but one real place [Bl Algorithm 4.1]. We find the algebra B = 
ramified only at vi and V2, generated by i, j subject to = w + 1, p = — 1, and 

For forms of parallel weight 2, Step 1 is trivial: we can take iiT = Q. 

Next, in Step 2 we compute ideal class representatives. The nontrivial class 
in Cl^(Zi?) is represented by the ideal b = {w^ — 2w — 6)Z^, which is principal 
but does not possess a totally positive generator, since sgn(— + 2w + 6) = 
(— 1, 1, —1) and there is no unit of l^p with this sign. We note that N(b) — 7. 

Next, we compute a maximal order O — 0(1); it is generated over Zj? by i 
and the element fc = (1 + (vj^ + l)i + Next, we find that the right O-ideal 

Jb generated by vo^ — 2w — 6 and the element (5 + {nP' +5)i + ij)/2 = 2 + 2i + k 
has nrd(Jt,) = b. 

Next, in Step 3 we compute presentations for the unit groups. We take the 
splitting 



B ^ M2(K) 




We find that F — r{l) is the free group on the generators a, /3, 71, . . . , 77 
subject to the relations 

7? = 72 = 7I = 74 = 7I = 76 = 77 = a/3a" V^Si • • • 77 = 1- 



w 



-1,-1 
F 



For example, we have 



2a = {w^ - 14) + (2^2 - 4w - 13)j + (-2w^ + 5w + 9)j + {-4w^ +8w + 26)ij. 



The groups F and Tf, have isomorphic presentations. In particular, we note 
that both r and Tf, have genus 1, so we conclude that dim 52(1) = 1 + 1 = 2. 

We illustrate the computation of Hccke operators with the primes p3 = {w + 
2)Zf of norm 3 and = {w + 3)Zp of norm 5. Note that p3 is nontrivial in 
CI'^{Zf) whereas p5 is trivial. 

Step Step 4 requires no work, since we work with forms of level (1). In Step 
5 we compute with a basis for cohomology, and here we see directly that 

H\r, Q) ^ Hom(r, Q) ^ Z/„ © Z/^ 

where fa,.fp are the characteristic functions for a and 13. We have a similar 
description for ff^(r[,,Q). 

Next, in Step 6 we compute representatives of the left ideal classes. For ps, 
for example, for /[i:o] C O we find that JbI[i:Q] = Ob{{w + 1) + i + ij) and for 
-^[1:1] C Of, we have Jb-^[i:ij = 0(ui + 1 — i + ij); we thereby find elements tTq, tt^ 
for a e P^(Fp3). For the generators 7 = a, ^ of and Ob, we compute the 
permutations 7* of P^(Fp3); we find for example that a* is the identity and 



[1:0] 



a 



0]''[l:0] 



with G Ob, namely. 



Ud[j^.Q^ = {7v? - 98) + {-2-iu? + 40w; + 167)i+ 

(-25w;^ + 59wj + 103),? + {-2w^ + Irw + 20)ij. 

We then write (5|^.q] as a word in the generators for of length 23. Repeating 
these steps (reducing a total of 64 units), we assemble the block matrix in Step 
7 as the matrix 

/O 2 0\ 
2 
2 
\0 2 0/ 

In a similar way, we find that Tp^ is the identity matrix. 

Finally, in Step 8 we compute the action of complex conjugation. Here we 
have simply = i (whereas /if, is more complicated), and thereby compute that 



W^\H = 



(I 1 \ 
0-100 
11 

\0 -1/ 



We verify that Wi 
/O 2 



commutes with Tp^ (and Tp^). We conclude that Tp^ \ 



2q) andTp, 



H+ 



1 
1 



We then diagonalize the space , which breaks up into two one-dimensional 
eigenforms / and 5, and compute several more Hecke operators: we list in Table 
1 below a generator for the prime p, its norm Np, and the Hecke eigenvalues 
ap(/) and ap(g) for the cusp forms f,g. 



p 


Np 






w + 2 




3 


2 


-2 


w + 2, 




5 


1 


1 


2 




8 


-5 


-5 


2io + 7 




9 


-2 


2 


w 




11 








— w — 


8 


17 


-5 


5 


TO — 3 




17 


-5 


-5 


2itP — 5w - 


- 10 


23 


2 


-2 


ttP — Sui - 


' 2 


25 


-9 


-9 


- 6 




29 


9 


-9 


li) + 4 




31 


-2 


-2 


2w^ - 3w - 


- 16 


37 


-3 


3 


— 2w - 


-9 


41 


-5 


5 


w'^ +w — 


3 


49 


-10 


10 



Table 1: Hecke eigenvalues for the Hilbert cusp forms for F = Q(iu) with — llui — 
11 = of level (1) and parallel weight 2 

We note that the primes generated by w and w — i are ramified in F . 

By work of Deligne the curves X = X{\) and X\, are defined over the 
strict class field of F, and Gal(F+/F) permutes them. We compute that 
F^ = F{\/ —ZuP- + 8w + 12). Therefore the Jacobian J/, corresponding to the 
cusp form /, is a modular elliptic curve over F+ with #J(Fp) = Np + l — a/(p) 
with everywhere good reduction. The form g is visibly a quadratic twist of / by 
the character corresponding to the extension F^ /F. 

Unfortunately, this curve does not have any apparent natural torsion struc- 
ture which would easily allow for its identification as an explicit curve given by 
a sequence of coefficients [H §4]. 

As a second and final example, we compute with a quaternion algebra defined 
over a quadratic field and therefore ramified at a finite prime. We take F = 
Q(y65), with Zf = + V65)/2]. The field F has #G\{F) ^ #C1+(F) = 2. 
We compute the space S — S'2(p5)''^""°" of Hilbert cuspidal new forms of parallel 
weight 2 and level ps, where p5 is the unique prime vu'Lf of norm 5. 

We compute that dim S = 10, and that the space S decomposes into Hecke- 
irreducible subspaces of dimensions 2,2,3,3. For example, the characteristic 
polynomial of Tp^ for p2 either prime above 2 factors as 

(T^ - 2T - 1){T^ + 2T~ 1){T^ + IIT'^ + 31T^ + 9). 

Remark 5. By the Jacquet-Langlands correspondence, the space S'2(p5)''^""''^ 
also occurs in the space of quaternionic modular forms for an Eichler order of 



level p5 in the definite quaternion algebra ramified at the the two real places 
of F and no finite place, and therefore is amenable to calculation by the work 
of Dembele and Donnelly. We use this overlap to duplicate their computations 
(as well as ours) and thereby give some compelling evidence that the results are 
correct since they are computed in entirely different ways. 
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